UNIQUENESS FOR KELLER-SEGEL-TYPE CHEMOTAXIS MODELS 
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Abstract. We prove uniqueness in the class of integrable and bounded nonnegative solutions 
in the energy sense to the Keller-Segel (KS) chemotaxis system. Our proof works for the fully 
parabolic KS model, it includes the classical parabolic-elliptic KS equation as a particular case, 
and it can be generalized to nonlinear diffusions in the particle density equation as long as the 
diffusion satisfies the classical McCann displacement convexity condition. The strategy uses Quasi- 
Lipschitz estimates for the chemoattractant equation and the above-the-tangent characterizations of 
displacement convexity. As a consequence, the displacement convexity of the free energy functional 
associated to the KS system is obtained from its evolution for bounded integrable initial data. 



I 

1. Introduction 
The classical Keller-Segel (KS) model for chemotaxis is the system 

dtn = kAu — xdiv (nVc), 



dtc = rjAc + On — 7c. 

Here, n is the number/mass density of a bacteria/cell population and c represents the concentration 
of a chemical attractant that can suffer chemical degradation and that it is produced by the cells 
themselves due to chemotactic interaction. The parameters k, x, V> ®-> 7 might be suitable functions, 
assumed to be constant in this simplified model. We can perform a time scaling and a suitable 
change of variables, that is r = Kt, p(x,r) = ^n(x,r/ k),v(x,t) = %c(x,t/k). The system is 
therefore reduced to 

( d t p = Ap- div(pVv), 
I edtv = Av + p — av, 

where a > and e > are constants (a = 7/77, e = k/w). In case e = 0, it restricts to the classical 
parabolic-elliptic Patlak-KS model 

J d t p = Ap - div(pVv), 
[ —Av + av = p. 

For e > 0, the natural free energy functional associated to the dynamics of the system (jl.ip is 

&e,a(P> v ) : = I (plogp-vp)dx + l [ (\Vv\ 2 + av 2 )dx. (1.3) 

In the case e = 0, corresponding to (|1.2|) . this Liapunov functional is at least formally equivalent 
to 

&oM : = (plogp- -vp)dx (1.4) 
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with the convention that v is obtained from the density p by v = & a ,d * P- Here, B a ^ denotes the 
Bessel kernel for a > or the Newtonian kernel for a = 0, for any dimension d. Therefore the role of 
the parameter e is to discriminate between parabolic-parabolic and parabolic-elliptic system. Note 
that the Liapunov functionals (jl.3p and (jl.4p are just formally equivalent since the L 2 -integrability 
of VB at d * p fails if d = 1, 2 and a = 0. Thus, even if the classical free energy writing and valid for 
all cases when e = is the one in (jl .4j) . we will prefer to work with the functional as in (jl.3p even 
if e = 0, with a suitable renormalization for the cases d = 1, 2 and a = discussed in Section 3. 

Our main objective is the uniqueness of certain solutions, for both systems (jl.ip and (jl.2p . Let 
us introduce the notion of solution for the Cauchy problems associated to (jl.ip and (jl.2p that we 
will consider in this work. We denote by ^(I ; m) the set of nonnegative densities over M. d with 
mass m and finite second moment, i.e., 

^# 2 (M d ;m) := \p G L 1 (R d ) : p > 0, / p(x) dx = m, \x\ 2 p(x) dx < +oo 

Definition 1.1. M^e say t/iai a weakly continuous map p G C w ([0, T]; ^(K ; m)) is a bounded 
solution to the Cauchy problem for (jl.2j) . uni/i initial datum p° G ^#2(K d ;ni) fl L°°(IR rf ) ; i/ 
ij pe L°°((0,T) x M d ) and H 2 /^) e L°°((0, T), L 1 (R d )), 

iij po = p° and the first equation of (jl.2p ZioWs in the sense of distributions on (0,T) x W 1 , 

where v t = B a ^ * pt for all t G [0, T], 
Hi) p t G W 1 ' 1 (M. d ) for C l -a.e. t G (0,T) and 

|Vp*(x)| 2 



o 



dxdt < +oo. (1.5) 



Definition 1.2. We say that the couple of functions (p,v), with p G C w ([0, T];^#2(R ;tn)) and u G 
L 2 ((0,T); W 1 ' 2 ^)), is a bounded solution to ([TT]) wi£a mi^a/ dainrn (p°,v°), p° G ^(I^; m) n 
L°°(M d ) and w° G W /1 ' 2 (IR d ), i/ 

I) p € L°°((0,T) x M d ) and |x| 2 pi(x) G L°°((0,T), L 1 (R d )), 

II) po = p° , the first equation of (jl.ip holds in the sense of distributions on (0, T) x M. d , and v is 
the unique solution to the Cauchy problem for the forced parabolic equation edtv — Av + av = 
p over (0, T) x M. d in the standard sense, with initial datum v°, 
III) the property Hi) of Definition II. ll holds. 

Let us emphasize that the main properties we need to get uniqueness of solution are the bound- 
edness of the densities and the Fisher information (II. 5p . They together imply that the velocity 
field of the continuity equation for the density p is a well defined object belonging to the right 
functional space, see section 2 for details. Moreover, the boundedness of the density implies that 
we have a uniform in bounded time intervals estimate on the quasi-Lipschitz constant of part of 
the velocity field. These are the basic properties that imply the uniqueness for bounded solutions. 
Let us finally mention that part of the strategy is related to the uniqueness of solutions to fluid and 
aggregation equations developed in [34], [27J, EHJ [HI [TJ [29] . The main novelty here is the interplay 
between the diffusive and the aggregation parts. The main results of this work are: 



Theorem 1.3. Let T > and let p u G ^(^nm) n L°°(R d ). Let pi,p2 be two bounded solutions 
on [0,T] x M. d to the Cauchy problem associated to (jl.2p . with initial datum p° . Then p\ = p2- 

Theorem 1.4. Let T > and let p° G ^2(M d ;m) n L°°(M. d ), v° G W^' 2 (M d ) n W 2 '°°(M. d ). Let 
(pi,v\) and (p2,V2) be two bounded solutions on [0, T] x M. d of the Cauchy problem associated to 
(jl.ip . with initial datum (p°,v°). Then (pi,v±) = (p2,V2)- 
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The proof of uniqueness as stated in Theorems 11.31 and 11.41 will be a consequence of a more 
general property: we will show that bounded solutions satisfies a strong gradient flow formulation 
by means of a family of evolution variational inequalities. This formulation is similar to the one 
for semi-convex functionals and implies a non-expansivity property of the distance between two 
solutions. This non-expansivity property yields uniqueness. All these results will be stated in 
Theorem 13.11 Moreover this formulation lead also to a relaxed convexity property of the energy 
functional as stated in Theorem 14.11 

There is a huge literature about the KS system and their variations, so we just restrict here 
to discuss the main results concerning about bounded solutions. In the classical parabolic-elliptic 
KS equation e = a = and d = 2, global in time bounded solutions in the subcritical case 
m < 8ir have been obtained joining the results in [12j [231 El- Actually, the global existence of 
weak solutions satisfying all properties in Definition 11.11 except the L°° bound was obtained in 
[12j while L°°-bounds in bounded time intervals can be obtained from the results in |24[ I14j . The 
same techniques could eventually be used to get local in time bounded solutions for all masses, 
although such a result is not present in the literature. Let us also mention the recent preprint [T7] 
in which the authors actually show that the L°°-norm of the solution decays in time like for the 
heat equation in the subcritical case m < 8ir for more restricted initial data. L°°-apriori estimates 
were obtained in the classical parabolic-elliptic KS equation e = with d > 2 and a > for small 
L d l 2 initial data in [201 [21] ■ These results together with similar arguments as in [12] to get the free 
energy dissipation property and thus the Fisher information bounds, could lead to the existence of 
bounded solutions in these cases. We emphasize that these L°° estimates show that the solution in 
bounded time intervals is bounded by a constant that depends only on the L°°-norm of the initial 
data, the initial free energy, and the final time. In particular, existence of bounded solutions is 
expected if p° £ L°°(M. d ), and this explains the presence of such an assumption in the previous 
definitions. 

Concerning the fully parabolic KS system, we find global in time solutions satisfying all properties 
stated in Definition 1 1 . 2 l except the L°° bounds in [15] for d = 2 and the subcritical mass case m < 8n. 
L°°-apriori estimates were obtained in [25] for the fully parabolic case but in bounded domains. 
It is reasonable to expect that this strategy should work for the whole space case, although it is 
not written as such in the literature. Results in higher dimensions concerning solutions with L°° 
estimates for small initial data can be found in [8] but estimates on the free energy dissipation are 
missing there. We finally refer to [111 \13\ [T4"l [22] for different results concerning the existence of 
solutions satisfying the boundedness of the Fisher information and/or the uniform bounds of the 
solutions for particular choices of e > 0, a > 0, and nonlinear diffusions. 

As mentioned before, Theorems 11.31 and 11.41 are based on the derivation of quasi-Lipschitz esti- 
mates for the chemoattractant v. This is the reason behind the additional assumption on the initial 
datum v° in Theorem 11.41 We will clarify the use of quasi-Lipschitz estimates of the chemoattrac- 
tant in the next section together with a quick summary of the main properties of optimal transport 
that we need in this work. Section [3] is devoted to show the main uniqueness results, derived from 
a more general property of bounded solutions for the Keller-Segel model. In fact, we will show 
that for bounded solutions we can obtain evolution variational inequalities. In Section [4] we show 
that these evolution variational inequalities lead to certain convexity of the associated free energy 
functional. In order not to break the flow of the argument, we postpone to Section [5] the rigorous 
derivation of the quasi-Lipschitz estimates of the elliptic and parabolic equations for v. In Section [5l 
we will also prove a strengthening of Theorem 11.41 with more general initial data. Finally, Section 
[6] is devoted to show how to adapt these arguments to Keller-Segel models with nonlinear diffusion. 
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2. Preliminary notions 



2.1. Some elliptic and parabolic regularity estimates. The proofs of our results are based 
on the technique used by Yudovich |34j for treating uniqueness in the case of incompressible Euler 
equations for fiuidodynamics. In particular, we exploit a quasi-Lipschitz property for the velocity 
field of the continuity equation for p in (jl.ip and (jl.2p . This property comes from the regularity 
that v gains being solution to the second equation in (|l.ip and (| 1 . 2 j) . 

Suppose first that v = Bq^*P- If p G L 1 nL°°(M. d ), by exploiting some estimates of the Newtonian 
potential, Vv satisfies the following log-Lipschitz property (see [6] and [301 Chapter 8], [32] and 
also [S]), 

\Vv(x) - Vv(y)\ < C\x - - log" \x - y\), 
where C is a suitable positive constant, depending only on ||p||x,i and ||p||l°° an d log - denotes the 
negative part of the natural logarithm function. As a consequence, we get the estimate 

\Vv(x)-Vv(y)\ 2 <C\{\x-y\ 2 ) (2.1) 
for some new positive constant C, where <p is the concave function on [0, oo) defined as 

/ \ ._) x\og 2 x if x < e' 1 '^, . . 

^ - \ x + 2(1 + V2)e~^ if x > e" 1 " A ^ 



Indeed, for large values of \x — y\ the estimate (|2.ip is quite obvious, since it is immediate to 
show that VBo.d * p is a bounded function in the whole space with a direct estimate using the fact 
that p G L 1 n L°°(M. d ). The log-Lipschitz property itself can be justified through standard elliptic 
regularity, as we will do in Section [5j 

Analogous facts hold if we consider the equation —Av + av = p, appearing in (jl.2p . or more 
general uniformly elliptic operators, so that we have 

Proposition 2.1. Suppose that p € L 1 r\L°°(M. d ) and a > 0. Then v = B a ( i*p satisfies the estimate 
(|2.ip . where C is a suitable positive constant, depending only on a,d, HpH^i^d), and \\p\\l°°{m.<i)- 

About the parabolic equation for v in (jl.ip . the quasi-Lipschitz property also carries over, since 
formally inequality (j2.ip translates in terms of the parabolic metric to 

\Vv(t,x) - Vv{s,y)\ 2 <C 2 <p({\x-y\ + \s-t\ 1/2 ) 2 ) Vx,yeM d , s,t € [0,T]. (2.3) 

Proposition 2.2. Suppose that p G L°°((0,T) x M d ) ; v° G Vy 2 >°°(]R d ) and a > 0. If v is the unique 
solution to the Cauchy problem for the parabolic equation dtv = Av — av + p (in the standard sense 
of convolution with fundamental solution), then v satisfies (j2.3p . where C is a suitable positive 
constant, depending only on a,d, \\v \\w 3 '°°(M. d )> an< ^ HpllL° c ((o,T)xR d ) ■ 

For a more complete insight into these properties, it will be convenient to recall some facts 
about the Zygmund class and its role in elliptic and parabolic regularity. However, in order not 
to introduce some not really necessary notation before the proof of our main results, we prefer to 
postpone the proof of Proposition [2TT1 and Proposition 12.21 to Section [5j Indeed, in Section [5] we will 
develop a more rigorous discussion about the log-Lipschitz estimates, and thanks to some refined 
parabolic regularity we will also prove a slight strengthening of Theorem 11.41 



2.2. Elementary notions of optimal transport. Given po,pi G ^#2(K , m), we define the 
Wasserstein distance between po and pi as 

W2(po,Pi)=( / \x - T(x)\ 2 p (x) dx 



where T is the unique optimal transport map between p^ and pi, that is, the map T : K — >■ 
K which minimizes J^ d \x — S(x)\ 2 po(x) dx among all the Borel maps S : M. d — > M. d satisfying 
5#Po = Pi- We recall that S#po = p\ means that J" Rd ip(x)pi{x) dx = J Rd ip(S(x))po(x) dx for every 
continuous and bounded function ip : M. d — > M. d . 

The Wasserstein geodesic between po and p\ is the curve s G [0, 1] i— > p s G ,y#2(R > ni) defined 
by the so-called displacement interpolation along the optimal transport map T between pq and pi, 
that is, p s := ((1 — s)i + sT)#po- I n particular, for any s, % := (1 — s)i + sT is the optimal map 
between po and p s and there holds W2(p r ,p s ) = \s — r\ W2(po, Pi)- 

We recall a formula for the differentiation of the squared Wasserstein distance along solutions of 
the continuity equation. Let t G [0, T] i— > pt G (K.^; tn) be a weakly continuous curve which is 
distributional solution of 

d t pt + div (£ t pt) = 0, 

for some Borel velocity field such that JcT H^* Hi 2 (iR d prffi 1 *) ^ < +oo. Then the curve is abso- 
lutely continuous with respect to the Wasserstein distance, [31 Theorem 8.3.1]. Then, for any 
p G ^ 2 (K d ;m), it holds 

1 A f 

--W 2 2 (p t ,p)= (Z t (x),x-Tt(x))p t (x)dx, for^-a.e. t€ (0,T), (2.4) 

2 at J R d 

where % is the optimal map between pt and p (see [3j Theorem 8.4.7, Remark 8.4.8]). 

Finally, let us recall an estimate relating the 2- Wasserstein distance and the H~ l norm proved 
in \Z7\ Proposition 2.8]. Given two nonnegative densities with the same mass pi,P2 G ^2(W i 'i m ) ^ 
L°°(R d ), there holds 

\\pi ~ Hl^-i^d) < max{||pi|| 00 , \\p2\\oo} 1/2 W 2 (pi,P2)- (2.5) 

Here we are letting if 1 (M rf ) be the space of Lebesgue measurable functions v : M. d — > M such that 
||Vw||^2(] R d) < +oo, so that H~ 1 (M d ) is defined by duality with functions having finite L 2 (M. d ) norm 
of the gradient only. By the way, we can also consider the space H 1 ^*) = W 1 ' 2 (R d ). In fact, from 
the proof in [271 Proposition 2.8] it is not difficult to see that the same estimate holds considering 
the i? _1 (R a! ) space given by duality with the full norm (|| Vv ||^2( K d) + ll^ll^^dp 1 ^ 2 - 

3. Bounded solutions as gradient flows: EVI and uniqueness 

The uniqueness Theorems 11.31 and 11.41 are consequences of a general result interpreting bounded 
solutions to (jl.ip (resp. (jl.2p ) as the trajectory of the gradient flow of the functional (|1.3|) (resp. 
(jl.4p ) in the appropriate metric setting. We prove that bounded solutions satisfy a family of 
evolution variational inequalities (EVI). Among different notions of gradient flow in metric sense, 
the EVI formulation is stronger than other formulations and typically corresponding to a convex 
structure, as in [3l Theorem 11.2.1] for the theory in the Wasserstein setting. 

Notation for the energy functional. Before giving the proof, we introduce some uniform 
notation for working with the full functional fjl.3|) even in the parabolic-elliptic case. Let p G 
^^(R^Siri) H L°°(M. d ). We are considering the free energy functional 

^£,a(p,v) := / (plogp - vp) dx + - / (|V-y| 2 + av 2 ) dx , 

defined for v being any W 1,2 (M. d ) function if s > 0. On the other hand, if s = it is understood 
that v is given by B a ^ * p. Therefore the parameter e only indicates if we are considering problem 
(jl.ip or (|1.2p . In particular, this writing of the functional as in (jl.31) is valid in general, even for 
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e = 0, except for two particular cases: e = a = and d = 1,2, as discussed in the introduction. 
In these two cases, we need to renormalize the free energy functional. Given p* G ^ 2 (W d ;m) a 
smooth and compactly supported density and v* = Bq^ * p*, we redefine fjl .3|) for e = a = and 
d = 1, 2 as 

&o,o(p,v):= [ [ p logp-v(p- p*)] dx+ \ f \X7(v-v*)\ 2 dx- \ [ p*v*dx. (3.1) 

Notice that V(u — v*) G L 2 (M. d ), as p — p* has zero mean, see [H[32] for more details. 

In the rest of this work, when referring to the free energy functional ^ e ,a (p> v ) > we will be using 
(jl.3p for any e > 0, a > 0, except for e = a = and d = 1,2 where the free energy functional is 
given by (|3.ip . 

Let us observe that now all the integrals involved, in the definition of a are well defined and 
finite for e > 0, a > and p, v as above. The negative part of the entropy term can be classically 
treated by the Carleman inequality, see for instance [9, Lemma 2.2] where the second moment 
bound on the density is used. The boundedness of the density controls the positive contribution of 
the entropy term together with the integrability of vp in case e > since v G W 1,2 (M. d ). For e = 
the integrability of vp in case a > is implied by the Newtonian potential case a = since the 
singularity of the Bessel potential at the origin is the same. The integrability for a = e = and 
d > 3 results directly from the Hardy-Littlewood-Sobolev inequality for the Newtonian potential. 
For a = e = and d = 1, 2 we use the behavior at infinity of the density p. Actually, a = e = 
and d = 1 is a trivial case since the Newtonian potential is given by Bqi{x) = \x\. For a = e = 
and d = 2 since log(e + |x| 2 )p G then vp G L 1 (M d ) using the logarithmic HLS inequality, 

see for instance [TO] . 

Notation for the ambient metric space. We let X £ := .^(B^; ni) X L 2 (M. d ) endowed with 
the distance 

D 2 (zi,z 2 ) = D 2 ((p 1 ,v 1 ),(p 2 ,v 2 )) = W$(pi,p 2 ) + e\\vi - v 2 \\ 2 L 2 



with the convention that Xq = ^^(R^jm) and Dq(zi,Z2) = W 2 {pi, p 2 )- Moreover, for z = p G 
Xq x L°°(]R rf ), ^(z) will be understood to be jFo,a(p>' u ) with v = Bq,^ * Pi as usual when e = 0. 

In the space X e the metric derivative of an absolutely continuous curve t i->- z% is denoted and 
defined by 

i /, / 1\ r Djzt+h^zt) 
\z \ D (t) = km , 

and it exists for £ 1 -a.e. i > 0. The local metric slope of the functional ^ £ ,a is defined by 

9^,ab W : = limsup '■ ' . 

These two abstractly defined objects are used to give the notion of curves of maximal slope in 
general metric setting, see [21 §3], [31 Chapter 1]. The main consequences of this gradient flow 
structure are summarized in the following result. 

Before stating the Theorem we define the function u : [0, +oo) —> [0, +oo) by 

oo(x) = y / mxip(m~ 1 x) 1 (3.2) 

where ip is defined in (I2.2j) . Moreover, given a fixed so > 0, we define a strictly monotone continuous 
function G : [0, +oo) — > [— oo,+oo) by G(s) := J" a * -^P)dr for s > and G(0) = — oo (we observe 

that G~ l : [— oo,+oo) — > [0, +oo) is surjective). 



Theorem 3.1. Lett \— > %i = (pt,Vt) be a bounded solution of problem (jl.ip fore > 0, starting from 
z° = (p°,v°) G X £ n (L°°(M d ) x (W 1 ' 2 (M d ) n Ty 2 '°°(M d ))) ; according to Definition lfe = 0, let 
zt = pt be a bounded solution to problem (jl.2p . starting from z° = p° G Xq n L°°(M. d ), according to 
Definition XY.W Then the three following properties hold: 

i) The evolution variational inequality (EVI) formulation: for any z = (p, v) G X £ C\ (L°°(W i ) X 
W 1 ' 2 ^)) (reduced to z = p G X n L°°(R d ) if e = 0), the map t (->• D 2 (z t , z) is absolutely 
continuous and there exists a constant C depending on ||p||Loo((o,T)xR d )> II/^IIl 00 ^) an d 
||t>°||^2,oo( R d- )J such that 

\^D 2 (z t ,z) < ^ a {z)-^ a {z t ) + Cuj(D 2 (z u z)) forC x -a.e. t G (0,T). (3.3) 

ii) The energy dissipation equality (EDE) in metric sense: the map t \- > ^ £ . a {zt) is locally 
Lipschitz continuous and 

j t &*M = -\\d,^ a \ 2 D {z t ) - \\z'\ 2 D (t) for£}-a.e. t G (0,T). (3.4) 

Hi) The following expansion control property: given another bounded solution t i— )■ Q, with 
initial datum £° in the same space of z° above, there exists a constant C , depending on 
ll/ 9 ll J L°°((o,T)xK d ) an d ll i;0 |lw /2 '° o (]R d ) (and the same quantities associated to (), such that there 
holds 

D 2 (z t ,(t) < G- 1 (G(D 2 (z°X°))+^Ct) for every t G [0,T). (3.5) 
Proof. We first introduce the auxiliary functional 

®e,a{p,v) := I (plogp-vp)dx, 

J]SL d 

for p and v being as in the definition of ^ £ia at the beginning of this section, so that 

^ £ , a {p, v) = $ £ , a {p, v) + - [ (\Vv\ 2 + av 2 ) dx 

and 

*o,o(p > u) = ^b ) o(p,«)- \ [ \V{v-v*)\ 2 dx + \ ( p*v*dx- [ p*vdx for d = 1,2. (3.6) 

^ JR d * JR d JR d 

The proof is organized in four steps. 

Stepl. Quasi- Lipschitz Estimate implies control of the evolution of the Wasserstein distance. - 
Thanks to the assumption (jl.5p , we learn that the Fisher information J Rd dx is finite for 

C l -&.e. t G (0,T). Let p G ^ (R rf ; m) n L°°(R d ). Exploiting the differentiability properties of the 
entropy functional, we can use the above-the-tangent formulation of displacement convexity to get 
for Z^-a.e. t G (0,T) 

p(x) log p(x) dx - / p t (x)\og p t (x) dx > / (Vp t (x),Tt(x) - x) dx, (3.7) 

JR d JR d 

where % denotes the optimal transport map between p t and p. We refer to [21 §3.3.1] for an 
intuitive proof of this fact, and to [31 Chapter 10] for the theory in full generality. In particular, 
the finiteness of the Fisher information of pt implies that the second term is finite, so that this 
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differentiation formula is meaning ful. If e > (resp. e = 0), let v € W 1,2 (K d ) (resp. u = B a4 * p). 
Take 



7 t := ® £}a (p,v) - $e,a(pt,V t ) + / (v(x) - V t (x)) p(x) dx . 

Using the notation xf := (1 — s)x + sTt(x), s £ [0, 1], and taking into account that 
v t (x)(p(x) - p t (x)) dx = (v t {T t (x)) - v t {x))p t {x) dx 

(v t (xl) - v t (x^))p t (x)dx = [ ^- [ v t (x s t )p t (x)dxds 
i ./n ds .had 



and (|3.7p . we obtain for £ 1 -a.e. t £ (0,T) 
It> / (Vpt{x),T t (x) - x) dx - v t (x)(p(x) - p t {x))dx 



(Vpt(x),T t (x) - x) dx -/ / (Vv t {x s t ),Tt(x) - x) p t {x)dxds 

1 JO JR d 

(Vp t {x) - pt(x)Vvt(x),Tt(x) - x) dx - / / (Vv t (xf) -Vv t (x),Tt(x) - x) p t (x)dxds. 
'■ Jo JR d 

Let us denote by TTf the last term in the right hand side above. The crucial point is to treat such 

term using the log-Lipschitz property of . Notice that, if e = 0, we are in the assumptions of 

Proposition 12.11 and we apply (|2.ip . where the constant C depends in principle only on (m, a, d 

and) the L°° norm of pt, which we are assuming to be uniformly bounded on (0, T). In the case 

e > 0, still by the uniform space-time L°° assumption on pt and the W 2, °° assumption on v°, we 

are in the framework of Proposition 12.21 so that we can apply the estimate (|2.3p . In this case the 

constant will depend also on (e and) ||v°|| H /2,oo( K d). Since tp is concave, we can also use the Jensen 

inequality, and letting p\ = x s t ^pt be the Wasserstein geodesic connecting p t and p we have 

\Ht\ <W 2 (p t ,p)£ n d \Vv t {x s t )-Vv t {x)\ 2 Pt (x)dx\ ' ds 

< CW 2 (p t , p) £ U v(\x s t - x\ 2 ) Pt {x) dx^j 1 ds (3-8) 

< V^CW 2 ( P t,p) J o ^(m^Wiipup^ds < V^CW 2 (pt,p)^(m- 1 Wi(p t ,p)) . 
The last inequality holds since geodesic interpolation ensures 

\x - x s t \ 2 Pt (x) dx = Wi( Pt , pi) = s 2 Wi(p t ,p) 



for all s S [0, 1] and since (p is non decreasing. We recall that the constant C in (|3.8p depends only 
on (e, a, d, the mass m and) the L°°((0,T) x M d ) norm of p and, in the case e > 0, the W 2,OQ (K d ) 
norm of v°. Inserting this in the estimate for It, we have for £ 1 -a.e. t £ (0,T) 

It> f (V P t(x)- p t (x)Vv t (x),Tt(x)-x)dx-Cu;(Wi(p t ,p)), (3.9) 

jR d 



where uj is the function defined in (|3.2p . Since pt satisfies the continuity equation 

d t pt + div (itPt) = with p£ t = -Vp t + pNv t 



and ([1.5p . the uniform L°° bound of pt implies that J Q T ||Ct||^2( R d pt .-^d\ dt < +00. Therefore 1 1-)- pt 
is absolutely continuous with respect to Wi and by (|2.4p 



1 d f 

o^W$(p t ,p) = / (Vpt(x) - pt(x)Vv t (x),T t (x) -x)dx for Z^-a.e. t G (0,T). 

2 dt J R d 

Inserting this into (|3.9p . and recalling the definition of It, we finally obtain 
1 d 



u Wi{ Pt ,p) < ®e,a{P>v) ~ * e , a {fit,Vt) + (V- V t )pdx + Cuj(W^p t ,p)) (3.10) 

z at j R d 
for ^-a.e. t G (0,T). 

Step EVI for the parabolic-parabolic case.- Recalling that v G W 1,2 (M. d ), observing that 
Av t G L 2 (R d ) for a.e.-i G (0, T) and using the elementary identity \a\ 2 — \ b\ 2 = \a — b\ 2 + 2(6, a — b) 
for every a, b G R fc , the variation of the second part of the functional (|1.3|) (that is, =^ eQ — $ ea ) 
can be written as 



If [\Vv\ 2 - \Vv t \ 2 + a(v 2 - v 2 )] dx 

2 JRd 



a , 



(av t - Av t )(v - v t ) dx + -||V(«t - v)\\ L2(Rd) + -\\v t - v\\ L2 



{p t - ed t v t ){v - v t ) dx + -\\V(v t - v)\\\ 2{Rd) + ^\\v t - v\\\ 2 

f— \ 7 S d — m2 lii, - — \ii2 ^ii — n2 

^ p t {v - v t ) dx + -— \\vt - v\\ L2m + -||V(u t - v)\\ L2(Rd) + -\\v t - v\\ L2 



Therefore, we deduce 

^e,a(p,V~) ~ ^e,a(Pt,V t ) = ® £ ,a(P,v) ~ ® £ ,a(pt,Vt) + I p t (v - V t ) dx 



[ Pt(i 



S d . . — 1 1 2 ^ii t 7 / — \ 1 1 2 ^ 1 1 — 1 1 2 



(3.11) 



(3.12) 



(3.13) 



Now, we use again (|3.10|) . leading to 

\j t W V t ~ ^Wl^) + ^^{pUp) ~ ^eAPiV) - ^e,a(pt,Vt) + C 'u (W| (ft , p)) 

By using the duality between H 1 and H , the Young inequality, and (|2.5p we have 

/" 11 

J {P ~ Pt){v ~ v t ) dx < \\p - pt\\ k „ l{Rd) \\v - v t \\ kl{w , d) < -\\p - Pt\\\- Hwd) + ^W^iv-vt) 

< ^QW 2 (p,p t ) + |||V(€; - v t )\\ 2 L2(Rd) , 

(3.14) 

where Q is the largest of the L°° norms of p and pt over the time interval (0,T). Taking into 
account that uj is given by (|3.2p and that \J mip(m^ 1 x 2 ) > x for every x > 0, combining (|3.13p and 
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II 2 



(|3.14p we get, up to introducing a new constant C, 

Idt"^ ~ ^'W d ) + \jt W ^ PU ^ ~ ^z^P^) ~ ^,a(Pt,V t ) + Cuj(W%(p t ,p)) ~ ^\\v t ~ vf L 2 

(3.15) 

for a.e. t G (0, T). The new constant C depends as usual on (e a, d, m and) ||/0||Lo°((o,T)xR d )) 
||« Q ||w r2 .° o 0R d )) 



Step 3: EVI for the parabolic- elliptic case.- When either d > 3 or a > 0, we can repeat the 
proof of the parabolic-parabolic case, letting e = therein and recalling that v is no more an 
arbitrary W 1,2 (M. d ) function but is given by convolution with p. In particular we arrive to the 
corresponding of (|3.13p . and the second line therein can now be estimated as follows. Using the 
inequality ||u||/r l -(]R < *) — llpll/f- 1 ^) f° r + av = p, a > 0, where the notation is |MI#i ( 
l|Vu||^ 2(Kd) + a|H|£ 2(R(i) (and using H 1 if a = 0), we get 

(v -v t )(p- p t )dx< \\v - v t \\ H i {Rd) \\p - p t \\ H ~i {Rd) < -\\v -v t \\ 2 H x iRd) + g IIP - Pt\\ 2 H , 
Moreover, recalling the estimate (|2.5p (which works both in H~ l and H~ l ) we have 

\\P- Pt\\H-\M,i) <Q W £(p>Pt)> 

for all t G [0,T], where Q is the largest of the L°° norms of p and pt over the time interval [0, T\. 
Inserting these estimates in (I3.13P we obtain 

\j t W i{puP) < ^AP)-^APt) + Cu>(WKpuP))i (3-16) 

for C l -&.e. t G (0, T), where the constant C depends only on e, a, d, m, ||/o||Loo((o,T)xR d )5 IIpIIl 00 ^)- 
In the case a = 0, d = 1, 2, we have to consider the functional in (|3.ip . By using the identity 

1 1 f 1 

2 ||V(tJ - v*)\\ 2 L2(Rd) - -\\V(v t - v*)\\ 2 L2(Rd) = J (pt - p*)(v - v t ) dx + -\\V(v t - v)\\ 2 L2{Rd) , 

with similar computations as in ()3.1ip . this time considering ^o,o(p> v ) ~ ^o,o(p 5 v ) as obtained from 
(|3.6p . we can still find (|3. 12[) and conclude obtaining again (|3.16p . 



Step 4-' Conclusion.- We are ready to prove the three points in the statement of the theorem. 
The proof of i) is a consequence of (|3. 15 j) for the case e > 0, and (|3.16|) for the case e = 0, taking 
into account that a > and that u(D 2 (zt, z)) > ui(W 2 (pt, p)) being to increasing. 

It is a standard fact that the gradient flow formulation in EVI sense implies the one in EDE 
sense in (|3.4p . Indeed, the proof of %%) follows from (13 .3p and (|3,5p and can be exactly carried out 
as in [2, Proposition 3.6]. 

The proof of (|3.5p still follows from (|3.3|) . Indeed we can apply [31 Lemma 4.3.4] (see also the 
argument of Theorem 11.1.4]) and obtain that for ^-a.e. t G (0, T) 



ld_ 
2~ds J 



1 d 



s=t 2 ds 



1 d 



D 2 {z s ,Cs) <^D 2 (z s ,(t) +-- r D\z t ,Q <2Cto(D 2 (z t ,Ct)). (3.17) 



8=1 



s=t 2 ds 

Here, C = max{Ci,C2}, where C% is the supremum on s G (0,T) of the constant in (|3.3p for zt 
with z = ( s , which is finite since £ G L°°((0,T) x M. d ), and C% is the same inverting z and £. The 
estimate (I3.17P implies 

^-D 2 (ztXt) < 4Cu(D 2 (z t ,( t )), for Z^-a.e. t G (0,T). 
at 
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Since the inequality 

y(t) <y(0) + 4C I u{y{s))ds 
Jo 

entails that y(t) < G _1 (G(y(0)) + 4Ct), we conclude. □ 

Proof of Theorems \1.3\ and \1.4\ The main theorems in the introduction are now a straightforward 
consequence of the expansion control Hi) in Theorem [37TJ Both Theorems follow from the inequality 
(133]) observing that G^C^O) + 4Ct) = G^f-oo) = 0. □ 



4. W-CONVEXITY OF THE FUNCTIONAL 

In this section we show another consequence of the EVI formulation of bounded solutions. For 
the functional ^ e ,a the following relaxed w-convexity along geodesies holds, see [18J for o;-convexity 
of functionals on measures. We assume that bounded solutions to (11. ip (resp. (11. 2p for e = 0) verify 
that for some T > 

\\pt\\L°°0i*) < Rt(p°,v°) := R(T,\\p \\ L oa m ,& £ , a (p ,v°f) , for ^-a.e. t G (0,T). (4.1) 

This assumption has been proved in several cases, see the introduction for more details. 

Theorem 4.1. Assume that bounded solutions for the evolutions (11. ip (resp. (II. 2D for e = 0) exist 
and verify gj). Then, for every z ^ 1 G X £ n (L°°(M d ) x (PF 1 ' 2 (R d ) n W 2 ' 00 ^))) (reduced to 
Xq D L°°(M d ) if e = 0) and every geodesic s G [0, 1] — > z s o/ i/ie space (X e , 13), connecting z° to z 1 
there holds for all s£ [0, 1] 

^.aM < (1 - s)^ e , Q (z°) + s^a^ 1 ) + fl T [(1 - s) ^(s 2 D V, z 1 )) + sw((l - s) 2 D 2 (z°, z 1 ))] , 
where Rt '■= max(i?y(z ), )). 

Proo/. We first remark that the set X £ n (L°° (M d ) x (VF 1 ' 2 (IR c( ) n W 2 ' 00 ^))) (resp. X n L°°(M d ) 
if e = 0) is geodesically convex. This is trivial for the part of the functional concerning v while for 
the density p we use the classical displacement convexity of all the L p norms |31j . Now, we take 
X £ n (L°°(lR d ) x (W l > 2 (R d ) n W 2 '°°(R d ))) (resp. X n L°°(E d ) if s = 0) as the set of initial data 
for the evolutions ([Lip (resp. (|1.2p ). 

Consider any z G X e n (L°°(IR d ) x W 1 ' 2 (M d )) (reduced to z G Xq x L°°(M d ) if e = 0). Taking 
into account that for a bounded solution t \— > D 2 (zt,z) is absolutely continuous and t i— >■ J^" ea (.Zt) 
is decreasing by ()3.4|) . from (|3.3p we obtain 

^D 2 (z t , f) - \d 2 {z , z) < t(& e , a (z) - & e , a (zt)) + C f io(D 2 {z r , z)) dr (4.2) 
1 1 Jo 

for all t G [0,T]. We denote by zf the bounded solution of (jl.ip or ()1 .2j) starting from the initial 
datum z s . We multiply by (1 — s) the inequality in (|4.2p for = zf and z = z° and we multiply 
by s the inequality in fj4.2j) for z^ = zf and z = z 1 . Summing up the two inequalities we obtain 

- s)D 2 (zl z°) + sD 2 (z s t , z 1 )) - - s)D\z s , z°) + *D V, z 1 )) 

< t ( (1 _ s ) J^ a (z°) + s^ a {z l ) - ^ a (zD) 

+ C((l-s) [ ' oo{D 2 {z s r ,z°))dr + s [ u(D 2 (z s r , z 1 )) dr). 
Jo Jo 
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Using the fact that s \— > z s is a geodesic, the right hand side is nonnegative, thus 

& e M) - (1 - s)^ £ , a (z°) - s^ £ta {z l ) 

<C((l-s)i f u{D 2 (z°,z°))dr + s- [ u{D 2 (z^z 1 ))dr). 
t Jo t Jo 

The lower semi continuity of 1 1— > and the continuity of r i— >■ D 2 (z%, z*), i = 0, 1 yield 

PeA*") ^ (! - «)-^,a(z°) + s-^aC* 1 ) + C((l - a)a;(Z> V, z )) + sa;(Z> 2 (z»> z 1 ))). 

Since s i-)- z s is a geodesic we have _D 2 (z s ,z°) = s 2 D 2 (z 1 , z°) and Z) 2 (z s ,z 1 ) = (1 — s) 2 D 2 (z°, z l ) 
and we conclude. □ 

5. A REFINED RESULT IN ZYGMUND SPACES 

This section is devoted to give a rigorous justification of the estimates stated in subsection 
2.1. We will also give a slight improvement of Theorem 13.11 and Theorem 11.41 by guaranteing a 
suitable quasi-Lipschitz estimate under a more general condition on the initial datum v . The right 
framework is that of Zygmund spaces. These classes of functions were introduced in [35], and they 
belong to the more general framework of Besov spaces. 

Zygmund estimates and log-Lipschitz regularity in the elliptic case. We first introduce the basic 
Zygmund class Ai(R d ) as the set of continuous bounded functions / over M rf such that 

\f( x )-2f((x + y)/2) + f(y)\ ^ 

SUp : : < +00. 

x,y€R d \ x ~ V\ 

It is well known that functions in the Zygmund class are in general not Lipschitz, possibly nowhere 
differentiable. Indeed, functions in Ai(R rf ) enjoy a log-Lipschitz modulus of continuity. Therefore, 
for any / £ Ai(IR rf ) there exists a positive constant C such that 



1/0*0 ~ f(y)\ < C\x - y\\\og\x - y\\ Vx,y G 



we refer for instance to [36J Chapter 2, §3]. More generally, following for instance [331 Chapter 5] 
we may define the class A n (R d ) for any n € N as follows. We let Ao = L°°(R d ) and we say that 
/ G A (R d ) belongs to A n (R d ), n > 2, if V/ G A n _i(R d ) or, equivalently, if / G W n - 1 '°°(R d ) and 
all the derivatives of / of order n — 1 belong to Ai(R rf ). In the usual notation of Besov spaces, A n 
corresponds to 00 . In this framework we have 

Proof of Proposition \2.1\ If a > 0, from the general theory on Bessel potentials (see for instance 
[331 Chapter 5, §3-6]) we learn that by convolution with the Bessel kernel B a ^ we indeed get two 
indices of regularity in A n spaces. Therefore, if p G L°°(R d ), we indeed get that v = B a ^ * p 
belongs to A2(R d ), and thus Vi> G Ai(R rf ) and ()2.ip follows. For the case a = we address to 
the references mentioned in Section [2] (it is also possible to directly check that Vi> G L°°(M' i ), and 
then the Newtonian potential behaves like the Bessel potential near the origin so that is also 
log-Lipschitz). □ 
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Zygmund estimates and log-Lipschitz regularity in the parabolic case. Let T > 0. In this section, 
let us denote Qt '■= (0, T) x M. d and then Qt := [0, T] x M. d . In the half d + 1 dimensional space, 
we consider the standard parabolic metric 



8((x,t), (y,s)) := max{|z - y\, y/\t - s\}. 

With respect to the parabolic metric, the definition of Zygmund spaces adapts as follows. We have 
Ao(Qt) := L°°(Qt), and Ai (Qt) is the space of continuous bounded functions / over Qt such that 
there hold 

\f(x,t)-2f((x + y)/2,t) + f(y,t)\ , \f(x, t) - 2f(x, (t + s)/2) + f(x, s)\ ^ 
sup : : + sup : < +00. 

t€[0,T] 0<s<t<T 

(5.1) 

Moreover, we say that / G L°°(Qt) belongs to A2(Qt) if 

gup \f(x,t)-2f(x,(t + S )/2) + f(x,s)\ < 

xeM d ,0<s<t<T I* ~~ s l 

and V/ G Ai(M d ). In particular, we see that / G A 2 (Qt) implies / G L oo ((0,T);W 1 ' oo (M. d )), with 
V/ satisfying (|5.1|) . so that finally / satisfies also (|2.3|) . 

When dealing with parabolic equations, it is suitable to consider spaces of functions defined with 
respect to the parabolic metric, since it is natural to deal with functions which have derivative up to 
order k with respect to time and 2k with respect to space. For classic results, we refer for instance 
to [IB] or to the monograph [2B] , where estimates are derived in Sobolev and Holder spaces of this 
kind, see Chapter 4 therein. 

In [16] we find that if the forcing term of the heat equation has bounded mean oscillation 
(BMO), still with respect to the parabolic metric, than the same holds true for second order space 
derivatives and first order time derivatives of the solution. This would be enough for deducing 
that first derivatives in space are in the Zygmund class with respect to the parabolic metric and 
that therefore they satisfy a log-Lipschitz estimate. The results in [16] deal only with null initial 
datum, but they can be generalized to more general data with suitable regularity requirements. 
Some extensions involving initial data in Zygmund classes are found in [H [23], based on direct 
estimates on fundamental solutions. Summing up, we have 

Proof of Proposition \2.tA Suppose that v is the solution (convolution with fundamental operator) 
of the forced heat equation dtv = Av + p. Suppose p G Ao(Qt) and v° G A2(lR d ). Then we have 
v G A2 (Qt)- See [16] for the case v° = 0, see [23|. Theorem 4] for a general result. If we consider 
the second equation of (11.11) with a > 0, the fundamental solution is just multiplied by a decaying 
exponential at infinity and the same result carries over. Therefore, a sufficient condition in order 
to have v satisfying fl2j| is v° G W 2 '°°(R d ), because W 2 '°°(R d ) C A 2 (R d ). □ 

This gives a rigorous justification of the assumptions on the initial datum of Theorem 11.41 
However a refined analysis shows that this assumption can be weakened, as we do next. 

Initial datum in A\(M. d ). We have to consider the weighted Zygmund space A 2 1 (Q T ), defined as 
the corresponding space A2(Qt), with the addition of a time weight which is divergent as t — > 0. In 
particular, locally in Qt functions in A^ 1 (Qt) have the same smoothness as the ones in A2(Qt)> but 
this regularity does no more extend to the closure of Qt- More precisely, by definition / G A 2 1 (Qt) 
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means that / G Ai(Qt), 

/7 \Vf(x,t)-2Vf((x + y)/2,t)+Vf(y,t)\ 
sup yt : : < +00 (5.2) 

x, v esi d \ x ~ y\ 

te[o,T] 

and the second finite differences of / and V/ with respect to time verify the corresponding estimates, 
as in the definition of ^{Qt), still with the addition of the weig ht tV2. 

Theorem 5.1. Let T > 0. Let p° G ^# 2 (R d ;m) n L°°(R d ) and v° G A x (M. d ) n W 1,2 (R d ). Let 
2^ = (p t ,vt) be a bounded solution on [0, T] x M d to £/ie Cauchy problem for (jl.ip . according to 
Definition 11.21 mt/i initial datum z° = (p°,v°). For any reference point z = (p,v) G (^(t ; m) D 
L°°(R a! )) x W 1 - 2 ^), t/ie genera/ EVI holds 

y^D 2 (z u z)< 3? £>a {z) - ^ £ ,a(z t ) + Ct- l ' 2 Lo(D 2 (z u z)) forC l -a.e. t€(0,T), (5.3) 

/or a constant C depending on ||p||£oo((o.T)xR d )> ll u °llA 1 (R d )j llpllL°°(K d )- 

Moreover the EDE (|3.4p holds, and the expansion control property holds in this form: given 
another bounded solution t H> £t as ak>i>e initial datum (° G (^(^im)!^! 00 ^)) x (Ai(M d )n 
W 1 ' 2 ^)) there is 



D\zt,(t) < G- l (G(D 2 (z°,(°)) + 8CVi) for every t G [0,T), (5.4) 

where C is a constant depending on ||/9||L°o(( ,T)xR d ) and \\v ||Ai(R d ) f ara ^ same quantities asso- 
ciated to Q). In particular, z = £ = £°. 

Proof. Since we are in the hypotheses of [23l Theorem 4], v belongs to (Qt), so that (15.2f) above 
holds for v and then, due to the log-Lipschitz regularity in the Zygmund class, we deduce 

\Vvt(x) — X7vt(y)\ < Kt~ 1 / 2 \x — y\\ log \x — y\\, (5.5) 

for all x G R d , t G (0,T), where K is a suitable constant depending only on T and the data. 
Notice that from the definition of A2 1 (Qt), it does not follow that Vt> G L°°(Qt)- However, in 
[231 Theorem 4] it is also shown that, still for v° G Ai(M d ) and p G L°°(Q T ), the solution v t of the 
parabolic equation also satisfies 

||V^(.)[|r»( R «i)<«'(l + |logi|), te(0,T), (5.6) 

again for some positive K depending on ||p[|i,oo(q t ) and 1 1 ^° 1 1 Ai • Taking (|5,6p into account, it is 
clear that (|5.5p can be improved (we do not relabel the constant) into 

\Vv t {x)-Vv t {y)\ < Kt-^lx-yKl + log- \x-y\). 

Thus we deduce the weighted analogous of (I2.3|) . that is 

C 2 

\Vv t (x) - Vv t (x)\ 2 < — <p(\x - y\ 2 ), (5.7) 

where C is a new suitable positive constant depending on the data and tp is defined in (|2.2p . 
Following the line of the proof Theorem 13.11 we reach the estimate (|3.8p for lit, which now has to 
be changed because we have to use (|5.7p . obtaining 

\Ih\ < Ct-^ 2 W 2 {p t ,p)^va V {m-^W 2 { Pu p)) = Ct~ l l 2 u{W 2 {p t ,p))- 

We can repeat all the other steps which lead to (|3.14p . obtaining the corresponding EVI with the 
additional weight t^ 1 / 2 , which directly lead to (|5.3p . We conclude as in Step 4 of the proof of 
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Theorem 13, It from (|5,3p . the EDE formulation (|3,4p follows, still referring to [2J Proposition 3.6]. 
Moreover, (j5.4f) follows by (|5.3p by 



d i D2 ( z utt) < ±Cr 1 l 2 u{D 2 {z t ,C,t)), for Z^-a.e. t G (0,T). 



Indeed the inequality 



y(t) < 2/(0) + 4C [ a-^uM*))** 
Jo 



implies that y(t) < G _1 ( G (y( )) + 8C\/F) as desired. Finally, the uniqueness result follows since 
G(0) = -oo and G-^-oo) = 0. □ 

6. The case of nonlinear diffusion 

We show next how to adapt our techniques to more general aggregation diffusion equations in a 
quite straightforward way. Let us consider the problem 

f dtp =div(pVP( / o))- div(pVv), 

< (6.1) 
I edfu = Av + p — av, 

to which we associate the functional 

%, a (p,v) := I {^{p)-vp)dx+\ I (\Vv\ 2 +av 2 )dx, (6.2) 

for all e > 0, a > 0, p G ^# 2 (K d ;m) n L°°(R d ), u 6 W 1 ' 2 (R d ), where := j?P{r)dr. We give 
the same restrictions as [31 §9.3], the first one being 

'l'(r) <i 
lim > — oo for some q > 



r^o ri d + 2 ' 

a property ensuring that L d VP(p) 7^ —00. Moreover, the crucial property to be satisfied by the 
new nonlinearity is the displacement convexity, that is the map r 1— > r d ^(r~ d ) is convex and non- 
decreasing on (0, +00). This notion, introduced in [31], is stronger than convexity and corresponds 
for C 2 functions to the inequality 

r -1 *(r) - *'(r) + r*"(r) > — r*"(r) V r G (0, +00). 

a — 1 

The more relevant cases correspond to nonlinear diffusion of power kind. Indeed, if 

m — 1 a 

the displacement convexity property holds. The case m > 1 (resp. m < 1) correspond to a slow 
diffusion (resp. fast diffusion) in the equation. On the other hand, the linear diffusion is recovered 
taking P(p) = log/3, it is seen that in this case functional (|6.2p is reduced, up to a constant, to 
(jl.3p . Finally, let us mention that the free-energy functional in the parabolic-elliptic case is similar 
to (jl.4p and given by 
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%, a (p,v) := / (*{p)--vp)dx, (6.3) 

JR d 1 

for p G ^2(1^;™) H L°°(R d ) and v = B a ^ * p. It can be written as (|6.2|) . taking into account the 
same renormalization as in (13. ip . to be done in the pathological cases e = a = and d = 1,2. 
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The notion of bounded solution is completely analogous to Definitions 11.11 and 11.21 both for the 
parabolic-elliptic and the parabolic-parabolic case. Indeed, the only point to adapt is the finiteness 
of the Fisher information, now rewritten into the generalized version 

f [ \VP(p t (x))\ 2 p t (x)dxdt< +00. (6.4) 

JO JR d 

Corollary 6.1. Theorem II .31 Theorem 1 1.4} Theorem 13.11 and Theorem 15. II hold for bounded solu- 
tions to (|6.ip . 

Proof. The displacement convexity property makes the internal energy functional p G ^2(R d ] m ) ^ 
L d ty(p(x)) dx convex along Wasserstein geodesies, as shown in [3, §9.3]. This in turn gives the 
possibility to write down a subdifferential inequality in Wasserstein sense (for a definition see [31 
§10.1.1]) as follows. Let p G Jfa(R d ;m) n L°°(R d ) be such that / Rd \VP(p)\ 2 pdx is finite. Then 

/ 9(p(x))dx- $(p(x))dx > / {VP(p(x)),T(x) -x)p(x)dx, (6.5) 

jR d JR d JR 

for any p G ^2(R d ',m), where T is the optimal transport map from p to p. Convexity and 
differentiability of functionals defined on probability densities, as the internal energy, are standard 
elements in the theory of Wasserstein gradient flows. For the proof of inequality (|6.5j) . which 
characterizes the vector V-P(p) as the Wasserstein subdifferential of the internal energy functional, 
we refer to §3.3.1] or to the general theory in |3j §10.4.3]. 

On the other hand, (16. 5ft can be used to generalize the proof of Theorem 13.11 Indeed, if (pt,vt) 
solves (16. 1|) according to our notion of solution, thanks to (16. 4b p t satisfies the identity (16. 5ft for 
almost any t. From this inequality, all the rest of the proof of Theorem 13.11 can be carried out. 
Indeed, with the same notation therein, we obtain the >C 1 -a.e. t G (0, T) inequality 

It ■= %,a{P, V) - %,a{pt,V t ) + / {v(x) - V t (x))p{x) dx > 

JR d 

f (p t (x)VP(p t (x))- Pt (x)Vv(x),Tt(x)-x)dx-Cu J (W 2 (p t ,p)), 

JR d 

for any p G ^# 2 (K d ; m)nL°°(IR d ) and any v G W l > 2 (R d ) if e > or v = B a4 *p if e = 0. This estimate 
substitutes (|3.9p in the proof of Theorem 13.11 The rest of the proofs is completely analogous. □ 
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